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Abstract 

We report a search for the charmless hadronic decay B° —> r/7r° with a data sample corresponding 
to an integrated luminosity of 694 fb -1 containing 753 x 10 6 BB pairs. The data were collected 
by the Belle experiment running on the T(4S I ) resonance at the KEKB e + e _ collider. We measure 
a branching fraction B{B° — > = (4.1+ 4 ' 5 lQy) x 10 -7 , where the first uncertainty is statistical 

and the second is systematic. Our measurement gives an upper limit of B(B° —> i]ir°) < 6.5 x 10 - ' 
at 90% confidence level. The signal has a significance of 3.0 standard deviations and constitutes 
the first evidence for this decay mode. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 
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Two-body charmless hadronic decays of B mesons are important for determining Stan¬ 
dard Model parameters and for detecting the presence of new physics [I]. The decay 
B° —* i)7T° proceeds mainly via a b —* u Cabibbo- and color-suppressed “tree” diagram, 
and via a b —» d “penguin” diagram [2], as shown in Fig. [lj The branching fraction can 
be used to constrain isospin-breaking effects on the value of sin 20 2 (sin 2 ct) measured in 
B —> tttt decays |3, Ej. It can also be used to constrain CP -violating parameters (C V , K 
and S V , K ) governing the time dependence of B° —> rj'K 0 decays [5j. The branching fraction 
is estimated using QCD factorization J 6 ], soft collinear effective field theory [7], and flavor 
SU(3) symmetry pj and is found to be in the range (2 — 12 ) x 10 '. 




FIG. 1. (a) Tree and (b) penguin diagram contributions to B° —> rj it 0 . 

Several experiments [9HT3], including Belle, have searched for this decay mode. The 
current most stringent limit on the branching fraction is B(B° —> ip r)° < 1.5 x 10 “ 6 at 90% 
confidence level (C.L.) |T3]. Here we update our previous result [12] using the full data set 
of the Belle experiment running on the T(4S') resonance at the KEKB asymmetric-energy 
e + e _ collider ra. This data set corresponds to 753 x 10 6 BB pairs, which is a factor of 5 
larger than that used previously. The analysis presented here also uses improved tracking, 
photon reconstruction, and continuum suppression algorithms. 

The Belle detector is a large-solid-angle magnetic spectrometer consisting of a silicon ver¬ 
tex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprising CsI(Tl) crystals. These de¬ 
tector components are located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside the coil (KLM) is instrumented to detect 
K® mesons and to identify muons. Two inner detector configurations were used: a 2.0 cm 
beampipe and a three-layer SVD were used for the first 123 fb -1 of data, while a 1.5 cm 
beampipe, a four-layer SVD, and a small-cell inner drift chamber were used for the remaining 
571 fb -1 of data. The detector is described in detail elsewhere [15, 46] . 

An ECL cluster not matched to any track is identified as a photon candidate. The timing 
of the energy deposited in the ECL must be consistent with the beam collision time identified 
at the trigger level. All photon candidates are required to have an energy greater than 50 
MeV. We reconstruct 7 r° —>• 77 decays by pairing together photon candidates and requiring 
that the 77 invariant mass be in the range 0.115-0.155 GeV/c 2 . This corresponds to ±3.5a 
around the nominal 7r° mass ms. Photon candidates in the end cap regions are required to 
have an energy greater than 100 MeV for 7 T° reconstruction. To improve the 7 r° momentum 
resolution, we perform a mass-constrained fit and require that the resulting y 2 be less than 
50. 

Candidate i) mesons are reconstructed via 77 —> 77 (r/ 77 ) and 7 —> 7 r + 7 r _ 7 r° (77 3 ^) decays. 
At least one of the photons in an 7/ 77 candidate must have an energy greater than 100 MeV. 
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To reduce combinatorial background from low-energy photons, r ^ 77 candidates are required to 
satisfy \Ei~E 2 \/{E\ +E 2 ) < 0.9, where E\ and E 2 are the two photon energies. Photons used 
for 77 77 reconstruction are not allowed to pair with any other photon to form a 7T° candidate. 
Candidate p 37T mesons are reconstructed by combining a 7 r° with a pair of oppositely charged 
tracks. These tracks are required to have a distance of closest approach with respect to the 
interaction point along the z axis (antiparallel to the e + beam) of \dz\ < 3.0 cm, and in the 
transverse plane of \dr\ < 0.3 cm. Pions are identified using information obtained from the 
CDC ( dE/dx ), the TOF, and the ACC. This information is combined to form a likelihood 
(£) for hadron identification. We require that charged tracks satisfy Ck/(C n + Ck) < 0.4, 
where Ck (/%) is the likelihood of the track being a kaon (pion). We reject tracks whose 
response in the ECL and KLM are consistent with that of an electron or muon. The pion 
identification efficiency is 91.2% and the probability for a kaon to be misidentihed as a pion 
is 5.4%. We require that the invariant mass of r / 77 and r/ 37r candidates be in the ranges 0.500- 
0.575 GeV/c 2 and 0.538-0.557 GeV/c 2 , respectively, corresponding to ±2.5cr and ± 3 . 0 cr in 
resolution around the nominal ij mass HU- For selected r) candidates, a mass-constrained 
fit is performed to improve the momentum resolution. 

Candidate B mesons are identified using the beam-energy-constrained mass, M bc = 
\]El e&ra - \p B \ 2 c 2 /c 2 , and the energy difference, A E = E B - E heam , where E heam is the 
beam energy, and E B and p B are the energy and momentum, respectively, of the B candi¬ 
date. All quantities are evaluated in the T(4S') center-of-mass (CM) frame. We require that 
events satisfy M bc > 5.24 GeV/c 2 and —0.30 GeV < A E < 0.25 GeV. We calculate signal 
yields in a smaller region M bc > 5.27 GeV/c 2 and — 0.21 GeV < A E < 0.15 GeV. 

Charmless hadronic decays suffer from large backgrounds arising from continuum e + e~ — y 
qq (q = u,d,s,c ) production. To suppress this background, we use a multivariate analyzer 
based on a neural network (NN) [18] . The NN uses the event topology and 5-flavor tagging 
information [19] to discriminate continuum events, which tend to be jetlike, from spherical 
BB events. The event shape variables include 16 modified Fox-Wolfram moments [20], the 
cosine of the angle between the z axis and the B flight direction, and the cosine of the angle 
between the thrust axis [21] of the B candidate and the thrust axis of the rest of the event. 
All quantities are evaluated in the T(45) CM frame. 

The NN technique requires a training procedure and, after training, achieves excellent 
separation between signal and background from the output variable Cnb ■ This variable 
ranges from —1 to + 1 : a value closer to —1 (+ 1 ) is more likely to identify a background 
(signal) event. The training samples used consist of Monte Carlo (MC) B° —> ipri 0 events for 
signal and an M bc -A5 sideband from data for background. The regions of 5.200 GeV/c 2 < 
M bc < 5.265 GeV/c 2 , A E < —0.23 GeV and A E > 0.17 GeV define the M hc -AE sideband. 
Independent samples are used to test the NN performance. The MC samples are obtained 
using EvtGen [22] for event generation and the Geant3 [53] package for modeling of 
the detector response. Final-state radiation is taken into account using the Photos [ 54] 
package. 

We require Cnb > — 0 . 1 , which rejects approximately 85% of continuum background 
events while retaining 90% of signal events. We subsequently translate Cnb to C' NB , defined 
as 


C'nb ~ l n 


( ri /^min \ 

%imax I 

°NB ~ °NBj 


(i) 


where C™ B = —0.1 and C“ff is the maximum value of Cnb obtained from a large sample of 
signal MC decays. This translation is advantageous as the C' NB distribution for both signal 
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and background is well described by a sum of Gaussian functions. 

After applying all selection criteria, 2% (7%) of events have more than one B° —» r^ 77 7T° 
(5° —>■ 7i 3n ir°) candidate. For these events, we retain the B° —> rj 7r° candidate with the 
smallest y 2 value resulting from the 7) or, if necessary, 7r° mass-constrained fits. According 
to MC simulations, this criterion chooses the correct B candidate 63% (77%) of the time for 
B° ->• ?7 77 7r° ( B° -)• r] 377 TT°). 

We calculate signal yields using an unbinned extended maximum likelihood fit to the 
variables M bc , AE, and C' NB . The likelihood function is defined as 


- Y Y i 


N 


C = e 


■n 


( 2 ) 


where N is the total number of events, AE\ C^ B ) is the probability density func¬ 

tion (PDF) of signal or background component j for event i, and j runs over all signal 
and background components. Y) is the yield of component j. The background components 
consist of e + e - —>■ qq continuum events, generic b —> c processes, and charmless rare pro¬ 
cesses. The latter two backgrounds are small compared to the qq continuum events and are 
studied using MC simulations. We find that no b —> c events pass our selection criteria. The 
charmless rare background, however, shows peaking structure in the M bc distribution, most 
of which arises from B + —y r]p + decays. 

Correlations among the fit variables are found to be small, and thus we factorize the 
PDFs as 

Vj(M hc , AE , C' NB ) = P,(M bc ) ■ Vj(AE) ■ V,{C’ NB ). (3) 

All PDFs for C' NB are modeled with the sum of two Gaussian functions. The M bc and 
AE PDFs for signal events are modeled with “crystal ball” (CB) functions [25] . The peak 
positions and resolutions in the signal M bc , AE, and C' NB are adjusted according to data- 
MC differences observed in a high statistics control sample of B° —> D°(— > K + n~7i°)'n' 0 
decays. This decay has four photons, as do signal decays, and its topology is identical 
to that of B° —* 7] 3n 7T°. The C' NB PDF of the continuum background is also adjusted by 
comparing data and continuum MC samples in the M bc sideband (5.200-5.265 GeV/c 2 ). The 
AE PDF for continuum background is modeled with a second-order polynomial, while the 
M bc PDF is modeled with an ARGUS function [9]. The M bc and AE PDFs for charmless 
rare background are modeled with one-dimensional nonparametric PDFs based on kernel 
estimation J2SJ• I n addition to the fitted yields Y-, the M bc and AE PDF parameters for 
continuum background are also floated, except for the end point of the ARGUS function. 
All other parameters are fixed to the corresponding MC values. To test the stability of the 
fitting procedure, numerous fits are performed to large ensembles of MC events. 

The signal yields obtained from the fits are listed in Table |T] The resulting branching 
fractions are calculated as 


B(B° ->■ r/7T°) 


Y ■ 

1 Slg 


yBB X 6 X 13ri 


(4) 


where U S i g is the fitted signal yield, N BB = (753±10) x 10 6 is the number of BB events, e is the 
signal efficiency as obtained from MC simulations, and B v is the branching fraction for rj —)■ 
77 or r) — y 0 H33- For the latter mode, e is corrected by a factor epio = 0.955 ± 0.015 

to account for a small difference in particle identification (PID) efficiencies between data and 
simulations. This correction is estimated from a sample of D* + —» D°(—> K~h + )'k + decays. 
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In Eq. Q, we assume equal production of B°B° and B + B pairs at the T(45) resonance. 
The combined branching fraction is determined by simultaneously fitting both B° —)■ ?7 77 7r 0 
and B° —> rj 3n n° samples for a common B(B° —y r\n r°). Projections of the simultaneous fit 
are shown in Fig. [2j 





FIG. 2 . Projections of the simultaneous fit: (a), (b) Mbc', (c), (d) A E; (e), (f) C' NB . Events 
plotted in Mb c (A E) are required to be in the signal region of A E (Mb c ), and C' NB >1.5 (0.5) for 
r 7 77 7 r° (? 737 rVr 0 ) decays. C' NB is plotted in the signal region of Mb c and A E. The left (right) column 
corresponds to 77 —> 77 (rj —> °) decays. Points with error bars are data; the (green) dashed, 

(red) dotted and (magenta) dot-dashed curves represent the signal, continuum and charmless rare 
backgrounds, respectively, and the (blue) solid curves represent the total PDF. 

The signal significance is calculated as 2 ln(£ 0 /£ max ), where £0 is the likelihood value 
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TABLE I. Fitted signal yield Uig, reconstruction efficiency e, r/ decay branching fraction B n , signal 
significance, and B° branching fraction B. The errors listed are statistical only. The significance 
includes both statistical and systematic uncertainties (see text). 


Mode 

Uig e(%) B v {%) Significance 0(10 ') 

B° —y ?7 77 7r 0 

SORIAS 18 - 4 39.41 

3.1 

r £+ 2.2 

5-6^2 .0 

B° %?+ 

0-5±i;® 14-2 22.92 

0.1 

0 . 2 l 2 | 

Combined 


3.0 

4 +L5 


when the signal yield is fixed to zero, and £ max is the likelihood value of the nominal fit. To 
include systematic uncertainties in the significance, we convolve the likelihood distribution 
with a Gaussian function whose width is set to the total systematic uncertainty that affects 
the signal yield. The resulting significance is 3.0 standard deviations; thus, our measurement 
constitutes the first evidence for this decay mode. A Bayesian upper limit on the branching 
fraction is obtained by integrating the likelihood function from zero to infinity; the value 
that corresponds to 90% of this total area is taken as the 90% C.L. upper limit. The result 
is B(B° —y ij-rr 0 ) < 6.5 x 10“'. 

The systematic uncertainty on the branching fraction has several contributions, as listed 
in Table [TT] The systematic uncertainty due to the fixed parameters in the PDF is esti¬ 
mated by varying them individually according to their statistical uncertainties. The result¬ 
ing changes in the branching fraction are added in quadrature and the result is taken as the 
systematic uncertainty. We evaluate in a similar manner the uncertainty due to errors in 
the calibration factors. The sum in quadrature of these two contributions constitutes the 
uncertainty due to PDF parametrization. We perform large ensemble tests in order to verify 
the stability of our fit model. We find a potential bias of —2.6%, which we attribute to our 
neglecting small correlations among the fitted observables. We assign a 3% systematic un¬ 
certainty for each reconstructed r; —» 77 or 7r° —>• 77 decay |2ZJ • The systematic uncertainty 
due to the track reconstruction efficiency is 0.35% per track, as determined from a study of 
partially reconstructed D* + —> D °(—* Kg7T + 7i~)Tc + decays. A 1.6% uncertainty (0.8% per 
pion) is assigned due to the PID criteria applied to charged pions in rj —» 7t + tt~tt 0 decays. 
We determine the systematic uncertainty due to the Cnb selection by applying different 
Cnb criteria and comparing the results with that of the nominal selection. The differences 
observed are assigned as the systematic uncertainty. The uncertainty due to the number 
of BB pairs is 1.3%, and the uncertainty on e due to MC statistics is 0.4%. The total 
systematic uncertainty is obtained by summing in quadrature all individual contributions. 

In order to check the reliability of the PDFs used for backgrounds, we fit the data in the 
Mbc sideband 5.24-5.26 GeV/c 2 , where the end point of the ARGUS function used for the 
continuum M^ c PDF is allowed to float. For all three distributions, Mbc A E, and C' NB , the 
MC PDFs give an excellent description of the data. We subsequently fit a sample of MC 
sideband events constructed with the same admixture of backgrounds as found in the data 
sideband and obtain signal yields consistent with zero. 

To check for potential nonresonant B° —>■ 777 r° and B° —y 7r + 7r~7r°n 0 decays, we relax 
the 77 mass requirement and plot the 77 and 7 r + 7 r _ 7 r° invariant mass distributions (Fig. [ 3 ]) 
for events in the M^ C -AE signal region. Significant peaks are observed for M 77 « M n and 








TABLE II. Systematic uncertainties on B(B° —> 777 T°). Those listed in the upper section are 
associated with fitting for the signal yields and are included in the signal significance and upper 
limit calculation. 


Source 

Uncertainty (%) 

PDF parametrization 

Fit bias 

+10.2 

- 9.2 

+0.0 

-2.6 

/r] —> 77 reconstruction 

6.0 

Tracking efficiency 

0.3 

PID efficiency 

0.6 

Cnb selection efficiency 

+2.1 

-2.2 

MC statistics 

0.4 

Nonresonant contributions 

+ 0.0 
-10.8 

77 ) 

0.5 

B(r) —> 7 r + 7 r _ 7 T°) 

1.2 

Number of BB pairs 

1.3 

Total 

+12.2 

- 15.9 


M 7r+7r _ 7r 0 ph M , as expected. The small sidebands indicate no significant contributions from 
nonresonant decays. We check this quantitatively by requiring that M 77 (M n + 7 r -„.o) be in 
the sideband 0.45-0.50 GeV/c 2 (0.56-0.58 GeV/c 2 ) and repeat the fitting procedure. We 
find 2.2l 3 ;| (—2.2l|; 5 ) signal decays, consistent with zero. To be conservative, we assign a 
systematic uncertainty due to B° —> 7771 ° nonresonant decays by appropriately scaling this 
fit result. 



FIG. 3. Distributions of (a) Af 77 and (b) M n+n _^ a invariant masses for events passing all selection 
requirements, except those for M 77 or M n+n _ v0 . 

In summary, we report a measurement of the branching fraction for B° —y rj 7 T° decays. 
We obtain 


B(B° T] 7 r°) = ( 4 . 17 J 7 ;?) 


x 10 


-7 
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where the first uncertainty is statistical and the second is systematic. This corresponds to 
a 90% C.L. upper limit of B(B° —>■ ijn 0 ) < 6.5 x 10~'. The significance of this result is 3.0 
standard deviations, and thus this measurement constitutes the first evidence for this decay. 
The measured branching fraction is in good agreement with theoretical expectations [5HH]- 
Inserting our measured value into Eq. (19) of Ref. [3] gives the result that the isospin¬ 
breaking correction to the weak phase </>2 measured in B —)■ tttt decays due to 7i 0 -r)-rj' 
mixing is less than 0.97° at 90% C.L. 
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